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Abstract 

The Y-chart is a powerful tool for understanding the relationship between various views (behav- 
ioral, structural, physical) of a system, at different levels of abstraction, from high-level, architec- 
ture and circuits, to low-level, devices and materials. We thus use the Y-chart adapted for graphene 
to guide the chapter and explore the relationship among the various views and levels of abstraction. 
We start with the innermost level, namely, the structural and chemical view. The edge chemistry 
of patterned graphene nanoribbons (GNR) lies intermediate between graphene and benzene, and 
the corresponding strain lifts the degeneracy that otherwise promotes metallicity in bulk graphene. 
At the same time, roughness at the edges washes out chiral signatures, making the nanoribbon 
width the principal arbiter of metallicity. The width-dependent conductivity allows the design of a 
monolithically patterned wide-narrow-wide all graphene interconnect-channel heterostructure. In 
a three-terminal incarnation, this geometry exhibits superior electrostatics, a correspondingly be- 
nign short-channel effect and a reduction in the contact Schottky barrier through covalent bonding. 
However, the small bandgaps make the devices transparent to band-to-band tunneling. Increasing 
the gap with width confinement (or other ways to break the sublattice symmetry) is projected to 
reduce the mobility even for very pure samples, through a fundamental asymptotic constraint on 
the bandstructure. An analogous trade-off, ultimately between error rate (reliability) and delay 
(switching speed) can be projected to persist for all graphitic derivatives. Proceeding thus to a 
higher level, a compact model is presented to capture the complex nanoribbon circuits, culminating 
in inverter characteristics, design metrics and layout diagrams. 
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The Gajski-Kuhn Y-chart (Fig. [T]) is a model which captures in a snapshot view, the 
essential considerations in designing semiconductor devices l| . The three domains of the 
Gajski-Kuhn Y-chart are on radial axes. Each of the domains can be divided into levels of 
abstraction, using concentric rings. At the top level (outer ring), we consider the architecture 
of the chip; at the lower levels (inner rings), we successively refine the design into finer 
detailed implementation: 

• Creating a structural description from a behavioral one is achieved through the pro- 
cesses of high-level synthesis or logical synthesis. 

• Creating a physical description from a structural one is achieved through layout syn- 
thesis. 
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FIG. 1. Gajski-Kuhn Y-chart adapted for graphene 



Building up such a multilevel view of graphene, a material of undoubted interest, is an 
ongoing process. The aim of this chapter is to touch upon just a few of the concentric circles 
to create such a viewing platform. 



18.1. THE INNERMOST CIRCLE THE ATOMISTIC VIEW 



Although there has been a lot of effort at this end from physicists and chemists, our 
aim here is more object-oriented - we wish to touch upon the GNRs, specifically, the edge 
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chemistry that leads to the observed lack of chirality and metallicity, then progressively 
outward along the Y-chart towards their two and then three terminal properties, and finally 
to an overall compact model describing its circuit level potential. 



18.1.1 Flatland: A romance in two dimensions 

Among elements on the periodic table, carbon is unique in that its sp^ planar and sp^ 
tetragonal bond energies are comparable, making it geometrically vary all the way from 3-d 
diamond to 0-d buckyballs 2|. Of this entire set, a candidate that combines impressive 
material properties (e.g. ultra-low effective mass), structural versatility, and amenability to 
industry-standard planar fabrication techniques, is undoubtedly graphene and its multiple 
derivatives, including carbon nanotubes (CNTs), bilayer graphene (BLG), epitaxial graphene 
(epi-G), strained graphene (sG) and graphene nanoribbons (GNR). 

The impressive electronic properties of CNTs and GNRs stem from their parent graphitic 
bandstructure . Without delving into the mathematics, it is important to rehash some of 
the salient features relating the bandstructure to the underlying chemistry. The hybridiza- 
tion of one s and two carbon p^. ,^-orbitals creates a planar honeycomb structure in a single 
graphene sheet, loosely resembling self-assembled benzene molecules minus the hydrogen 
atoms. The crystal structure can be described as a triangular network with a two-atom 
dimer basis, whose vr electrons hybridize to create bonding-antibonding pairs that delocalize 
over the entire crystal to generate conduction and valence bands. However, since the two 
basis atoms and the orbitals involved are identical, we get a zero-band gap metal with a 
dispersion resembling photons, albeit with a much lower speed. The resulting low energy 
linear dispersion corresponds to a constant slope and thus a constant velocity v = dE/hdk. 

The unique bandstructure of graphene contributes to its amazing electronic properties . 
Because the Fermi velocity is energy-independent, the cyclotron effective mass of graphene 
electrons, m* = hkp/vF is vanishingly small at low energy [kf — 0, vp being constant at 
roughly 10^ cm/s). Furthermore, the two bands are derived out of symmetric and antisym- 
metric (bonding and antibonding) combinations of the two identical dimer atoms, creating 
a two component pseudo-spinor out of the two Bloch wavevectors, with their ratio being 
just a phase factor e*^, where tan^ = ky/k^ relates electron quasi-momentum components in 
the graphene x-y plane. The reversal of phase between the forward and backward velocity 



vectors suppresses ID acoustic phonon back-scattering, allowing only Umklapp processes in 
confined graphitic structures such as CNTs and GNRs. The combination of low mass m* 
and high mean free path A ultimately leads to very high mobilities fi = qX/m*VF, with a 
record room-temperature value at 230,000 cm'^/Vs {4] for suspended graphene sheets. 

In the following sections, we discuss two bandstructure related issues that arise when we 
attempt to pattern or modify graphene to generate gapped or confined planar structures - 



the absence of metallicity and chirality dependent bandgaps in multiple experiments 
(sections 18.1.2-18.1.4), and 

the increase in effective mass as a gap is progressively opened, arising from fundamental 
asymptotic constraints (section 18.1.5) that are expected to persist even for very pure 
samples. 



18.1.2 Whither metallicity? 

The bandstructure of CNTs is by now, a textbook homework problem. Since the con- 
duction and valence bands of its parent graphene structure touch precisely at the vertices 
of its hexagonal Brillouin zone, the Fermi wavelength of undoped graphene corresponds to 
a unique electron wavelength Xp = 3\/3R/2, where R is the C-C bond distance. The impo- 
sition of periodic boundary conditions along the circumference of a CNT filters out many 
allowed electron wavelengths and allows only modes that have integer ratios of irD/X [D 
being the tube diameter), so that a particular chirality (i.e., a wrapping topology) may or 
may not support Xp needed to sustain metallicity. Accordingly, one can derive selection rules 
on paper, or using simple 1-orbital orthogonal nearest neighbor tight binding, that matches 
experiments quite well. In any random array of CNTs, roughly a third are expected to be 
metallic, and impressive progress has been achieved in sorting them out from their semi- 
conducting counterparts. Even in extremely narrow CNTs with strong curvature-induced 
out-of-plane hybridizations, the anomalous bandgaps are well captured by non-orthogonal 
tight-binding formulations such as Extended Hiickel Theory [sl. The bandgaps of CNTs 
bear relatively few surprises, at the end of the day. 

Life is more complicated when dealing with GNRs. Indeed, it seems reasonable to expect a 
chirality dependence to arise for GNRs, simply replacing the periodic circumferential bound- 
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:^IG. 2. (Left) Tight-binding 1-orbital calculations show three chiral curves, one of which is metallic 
^. In contrast, data from (center) Philip Kim's group 7] and (right) Hongjie Dai's group ^ 
show a single chirality free curve with no metallic signatures. 

ary conditions across CNTs with hard-wall boundary conditions across the GNR width. A 
few details may change, for instance, we now fit half-wavelengths rather than whole wave- 
lengths, and the edges are not completely opaque to electrons tunneling outwards so that 
the boundary conditions are more 'diffuse'. The quantization condition will roughly corre- 
spond now to integer ratios for (W + R\^)/\ (accounting for two unit cells outside for the 
wavefunction to vanish). However, one-orbital tight binding would still predict three chiral 
classes in GNRs, one of which is strictly metallic as in CNTs. Experimentally however, no 
chirality dependence is observed for GNRs, nor are any GNRs observed to be strictly metal- 
lic at low temperature (Fig. [2]). Regardless of wrapping vector, GNRs wider than 10 nm 
are quasi-metallic while narrow ones semiconducting. We thus have a startling disconnect 
between simple theories and experiments 



18.1.3. Edge chemistry — benzene or graphene? 

The disconnect between naive expectations and observations arises from the boundaries, 
which ultimately impact the quantization rules behind patterned GNRs. Specifically, we ar- 
gue that edge strain and roughness are the main factors behind the disconnect. Doing justice 
to such effects will need a proper bandstructure that can capture atomistic chemistry and 
distortion. While Density Functional Theory (DFT) within the LDA-GGA approximation 
captures these effects well, we used DFT primarily for structure evaluation, and resorted to 
Extended Hiickel theory (EHT) fitted to bulk graphene to explore the low-energy bandstruc- 
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FIG. 3. C-C bond length comparisons show that relaxed armchair GNR edges lie between benzene 
and double bonds, enjoying only a partial resonant hybridization. 

tures. We relaxed the hydrogenated edges of armchair GNRs using LDA-GGA and found 
a bonding environment distinct from bulk graphene. While the inner carbon atoms have a 
bond-length of 1.42A, the edges tend to dimerize and see a 3.5% strain associated 

with a reduction in bond length to 1.37A. Thereafter, we employ non-orthogonal basis sets 
in EHT to capture the effect of the edge chemistry on the low-energy electronic structure. 

There is an appealingly simple explanation for the observed bonding chemistry. Since the 
edge atoms are connected to hydrogen on one side and carbon on the other, the difference in 
electronegativity tends to move the edge C-C atoms closer to a benzene structure (Fig. [3]). 
However, the unequal bonding environment at the edge disallows any resonant hybridization 
that evens out the double-single bond distribution in aromatic rings, so that the edge rings 
in GNRs break into 'domains' with nearly intact double bonds at the edges and slightly 
expanded single bonds towards the bulk end. Since benzene is semiconducting, the 3.5% 
strain at the dimerized edges increases intradimer overlap but reduces interdimer overlap, 
effectively opening a bandgap by 5%. The obvious consequence is that all armchair GNRs 
become strictly semiconducting, in sharp contrast to their CNT counterparts. In contrast, the 
edges of zigzag GNRs have a lateral symmetry that makes them resistant to dimerization. 
In fact, the inward motion of the C-C edges away from the hydrogens shrinks both bond 
lengths equally, making zigzag edges more conducting. 

To calculate the impact of the relaxed bonds on the electronic structure, we calculated the 
density of states of a uniformly wide armchair GNR with edge relaxation using Extended 
Hiickel Theory. While EHT has been used widely to study molecular properties, it has 
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FIG. 4. (a) Open carbon bonds at the edges introduce edge-states (shaded) in the DOS. Spatial 
resolution of those eigenstates around the Fermi energy confirms the electron wavefunction localized 
at the armchair edges, (b) When open carbon bonds are hydrogen terminated, those edge-states are 
removed, (c) Applying EHT to GNR dispersion relation across a range of sub-lOnm armchair edge 
widths finds an oscillating bandgap. (d) Relaxation of edge bonds that are hydrogen terminated 
widens the energy bandgap for 3p and reduces the gap for 3p + 1 GNRs. Eg vs width results are 
within range of experimental data points 8|] and also in agreement with DFT predictions. 

also been extended to describe bulk semiconducting bandstructures using localized Wannier 
like non-orthogonal basis sets that still retain their individual orbital properties. Through 
extensive tests on both graphene and silicon, we have found that EHT accurately captures 
both bulk bandstructures as well as surface and edge distorted bandstructures The 
result of our simulation is shown in Fig. |H The role of edge passivation is shown at the 
top, where we can see the explicit removal of edge induced midgap states by hydrogenation. 
The bottom panels show the role of edge strain. In contrast to pz based nearest neighbor 
one orbital tight-binding theory, a small bandgap opens jsl. While CNTs have precise 
periodic boundary conditions along their circumference, the edge atoms do not provide an 
exact hard wall boundary condition, as the electrons tend to tunnel out into the surrounding 
region. In the presence of edge relaxation, the gap increases because of the aforementioned 
dimerization, removing any semblance of metallicity from the bandgap vs width plots. 

While EHT explains the removal of metallicity, compact models prefer a suitably cal- 



ibrated orthogonal tight-binding model, with the edge chemistry through beyond nearest 
neighbor interactions. It is not clear if this reproduces the Bloch wavef unctions, but they 
do seem to capture the overall density of states. Fig. shows a comparison between a (9,0) 
EHT DOS and a (9,0) tight-binding Density of States (DOS) (parametrized independently 
by CT White) We will use this formulation for our simpler compact models described 

later. 
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FIG. 5. Atomistic EHT plot (left) of the DOS of a (9,0) relaxed armchair GNR, compared with 
a next-nearest neighbor tight-binding fit (right) from C. T. White at NRL. From the pz sector of 
the EHT Hamiltonian, we estimate the nearest neighbor coupling, third nearest neighbor coupling 
and the shift in coupHng at the 3.5% strained edges as ti = —3.31 eV, = —0.106 eV and 
Ati = —0.38 eV. These extracted parameters compare quite favorably with C. T. White's tight 
binding parameters, ti = —3.2 eV, = —0.3 eV and Ati = —0.2 eV. [9| 



18.1.4. Whither chirality? 

While we can explain away the lack of metallicity through the preponderance of edge 
strain, why do we not see the three chiral curves in a plot of experimentally measured 
bandgaps versus ribbon widths? The primary reason, we believe, is roughness at the edges, 
which tends to wash away such chiral signatures. Currently, line edge roughness along a 
GNR edge is an unavoidable consequence of lithographic and chemical fabrication techniques 
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12 1 13|. Unzipping CNTs via ion bombardment produces the smoothest GNR edges to date 
14j |. However edge fluctuations even on the scale of a single atom can degrade transmission 
probability of modes near band edges, which have implications we expound in later sections 
on material and device characteristics such as mobility, subthreshold swing, and ON-OFF 
ratio. 
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FIG. 6. (Top) Two kinds of roughness include variation in GNR width and a width dislocation 
across a slip line, maintaining the same width (Bottom) Transmission plots show that for either 
case, segments with the larger bandgap filter out the rest of the segments, thereby promoting the 
chiral curve with the highest bandgap, in agreement with experiments and EHT calculations [lo| . 



We have studied a wide spectrum of line edge roughnesses that can ultimately be classified 
as either a width modulation or width dislocation. Modulation in width along the armchair 
edge has a corresponding modulation in bandgap that follows an oscillating inverse square 
law relation between bandgap and width, spanning the three chiral curves. Meanwhile a 
width dislocation is an in-plane displacement in the GNR that energetically sees the same 
bandgap at interface of the slip dislocation, albeit with localized interfacial states. We 
introduced edge roughness in our geometries using a Gaussian white noise that adds or 
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removes an integer number of dimer pairs with a correlation length 111- The results of 
our EHT simulations with a statistics of rough edges explain why measured data are seen 
to cluster around one the 3p + 1 chirality curves. 

An E-k relation cannot be rigorously defined for a structure without periodicity, so we 
focus instead on the transmission bandgap of the entire structure, calculated using the 
Non-Equilibrium Green's Function (NEGF) formalism in its simplest, Landauer level im- 
plementation. The plotted transmission of the rough segment (Fig. [6]), sandwiched between 
two bulk metallic contacts, is shown for a random mixture of (7,0) and (9,0) chiralities 10 |. 
The rough GNR can be viewed as a random mixture of individual GNR segments with well 
defined chiralities and bandg long as the correlation length is larger than the 

electron wavelengths. In the large majority of our simulations, we find that the transmission 
bandgap for the rough GNR follows the largest bandgap of the individual segments (in this 
case, the (7,0) segment). The larger bandgapped regions effectively filter out the incoming 
currents from the rest of the crystal, explaining the observed predilection towards the 3p -|- 1 
chiralities. There are a few outliers however, corresponding to accidental situations where 
the segments with large bandgap are thin enough to become transparent tunneling. As we 
increase the frequency of roughness by bringing the correlation length below the electron 
wavelength, we expect to see more tunneling events. Such high frequency roughnesses, as 
well as inelastic scattering, will mix chiralities more thoroughly and are predicted to spread 
out the data more evenly over the three chiral curves rather than clustering them around 
the dominant (3p-|-l) segment. 

The upshot of this analysis is that for reasonably smooth GNRs with Lm > A^, we 
expect the GNR width to be the sole arbiter of metallicity, so that wide GNRs are metallic 
while narrow (< 10 nm) GNRs are semiconducting |8|. While the analyses leading to 
this conclusion assumed some simplifications (e.g. ignoring structural modifications due to 
back-bonding with substrates for instance), experiments seem to support this conclusion. 



18.2. THE NEXT CIRCLE TWO TERMINAL MOBILITIES AND I-VS 

Our next circle would move on from the material parameters to electronic properties 
such as its carrier density, mobility, conductivity, and ultimately its current-voltage (I-V) 
characteristics. 
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18.2.1. Current-voltage characteristics (I-Vs) 



J T{E)M{E)[f,{E)-f2{E)]dE (1) 



The Landauer expression gives us a convenient starting point for the current through any 
material, 

-I 

where T ^ A^c/ (A^c + L) is the quantum mechanical transmission per mode, that relates its 
scattering length Xsc with its length L. The number of modes M = '-feffD{E), D{E) being 
the density of states, and the effective injection rate is given by I /'Jeff = I/71 + I/72 + 
^/ich, where 71^2 are the broadenings from the contacts, and •jch = hv{E)/L is the intrinsic 
transport rate in the channel. Assuming the contact broadenings are large so that the rate 
limiting step is 7c/i, we can replace 7e// ~ 7c/i. 

The band dispersion of graphitic materials, ranging from epi-G to sG, BLG, CNTs and 



GNRs are all described by a universal formula ISj 



E = ±^Ely + mvlk^ (2) 

where the band-edges are at Ecy while the high energy velocity in the linear regime is 
~ 10^ cm/s. From the dispersion, we can readily extract the 2D density of states and 
band velocities 



/ 2WL\ 



D{E)=\—-,\\E 



e{E - Ec) + e{-Ev - E) 



v{E)=vo^l-Ely/E^ (3) 

There is an additional energy dependence in the scattering length A^c- For ballistic channels, 
this is energy-independent, while for charge impurity and edge roughness scattering, A^c oc E, 
while for acoustic phonon scattering, Xgc oc 1/E. The actual dependences are a bit more 
complicated, but these are reasonable approximations to adopt. 

The algebra becomes particularly simple if we ignore the energy- dependence of Xsc- We 
can then do the Landauer integral, leading to 



Sq XW 



h \ nhvoL 



(4) 



where the shape function Iq depends on the current flow regime. Assuming we start with 
an n-doped graphene with a bandgap, we get 
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where /xi = -E^ and /i2 = -Ef — Q'^d- The expressions can be further simphfied. In the 
hnear regime, the current looks hke 

2 
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where Go = Q'^/h, the number of modes M ^ 2W/{Xf/2), and the Fermi velocity vp 
foA/l — El/Ep. The saturation current 

Ep 
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(7) 



while the band-to-band tunneling current at high bias varies quadratically as 

qVD 



W~4GoM| — ]Vd 




2Ef 



(8) 



Fig. [7] shows typical I-Vs based on the above formula. These results agree with more 
involved, atomistic models for EHT coupled with non-equilibrium Green's function based 
simulations 16|. The current shows a point of inflection at the Dirac point, which is shifted 



by the gate bias (bandgaps would give more extended saturating regions, as we will see for 
our three terminal I-Vs later on). The subsequent rise in current is indicative of band-to- 
band tunneling. Furthermore, a prominent I-V asymmetry, consistent with experiments on 
SiC, can be engineered into our I-Vs (Fig. |H]) readily by shifting the Fermi energy to simulate 
a charge transfer 'doping [l^]' of 470meV through substrate impurities, back-bonding and/or 
charge puddle formation with SiC substrates. A mean-free-path(Asc) that varied inversely 
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FIG. 7. A typical I-V output shows how the I-V tends to saturate at the Dirac point even without 
a bandgap. The shift in the Dirac point indicates the Laplace potential drop along the channel, 
eventually leading to band-to-band tunneling. 

with gate voltage was implemented in the left figure in Fig. |H1 For n-type conduction A^c 
ranged between 18nm to 40nm and 20nm to 31nm for p-type. Typically we would expect 
at least lOOnm for low bias conductance and down to lOnm as the biasing approaches the 
saturation and band-to-band tunneling regions. Chosen A^c represent an average scattering 
length for the different regions. A more accurate model for scattering is necessary and will 
be developed in future works. In contrast, Si02 seems to dope the sheets minimally and 
the measured I-Vs show the expected sjTiimetry between the electron and hole conducting 
sectors. 

From the low-bias I-Vs, we can now extract the conductance G = Iiinear/VD, thence the 
sheet conductance ag using G = agW/L, and finally the mobility /i^ using cr^ = qUgfis, where 
Ug is the sheet charge density related to the Fermi wavevector as kp = i/vrnj. Let us focus on 
the mobility first, keeping in mind that the effective mass m* in GNRs is energy- dependent. 

18.2.2. Low bias mobility-bandgap tradeoffs: asymptotic band constraints 

Graphene's linear dispersion is known for contributing to an ultra-high mobility, but of- 
ten overlooked is its origin in the low bandgap that ultimately hampers its ON-OFF ratio 
as an electronic switch. The carrier velocity {v = 1/hdE/dk) fundamentally saturates to 
Vq = 3aot/2h ~ lO^cm/s, which forces its high energy bandstructure to a linear form regard- 
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FIG. 8. (Left) Theoretical and (Right) experimental I-Vs for graphene. The calculations on the 
left assume a 'doping' of the sheet by a charge density that shifts the i^-point relative to neutrality. 
We also assume an inverse relation between the scattering length A^c and the applied voltage on 
the n-side, consistent with scattering by charge puddles associated with the above doping charge. 
The data on the right are for graphene on SiC, where charge puddles and/or back-bonding are 
expected to transfer a net charge density to the sheet {l8|. 

less of bandgap size. Regardless of the mechanism of bandgap opening, or the particular 
progeny of graphene that we are looking at (epi-G, sG, CNT, GNR or BLG), the bandstruc- 
tures are always writeable as E{k) ~ ±^y [Eg/^Y + {hvoky {2]. Such an intimate relation 
between bandgap and dispersion connecting ultimately to its conduct ion/ valence band effec- 
tive masses (as opposed to mid-gap tunneling effective mass) is unique in materials science. 

An extended bandgap constrained by the high energy velocity saturation localizes carriers, 
which shows up as a decrease in curvature at the band edges. It is easy to show from the 
above dispersion that the effective mass at each band-edge satisfies 



rrin 



(9) 



indicating that the kinetic energy gained by the electrons and holes upon bandgap opening 
is taken ultimately from the corresponding crystal potential (Fig. [9]). 

For transport considerations corresponding to high bias electrons and holes, we need to 
generalize the concept of effective mass to points away from the band-bottom, using m* = 
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FIG. 9. Scaling of graphite effective masses shows that increasing the bandgap increases the mass 
m* = pf/vf due to the decrease in average curvature arising from a pinning of the E-k at high 



energy values 



p/v = hk/[l/h{dE/dk)]. This expression reduces to the usual dependence on curvature 
near the band-bottoms upon using L'Hospital's rule with /c — )■ 0. In other words, effective 
mass and carrier velocity must be treated as energy dependent variables instead of material 
specific constants. From here, we can then extract the mobility fi = qXsc/p, where carrier 
momentum p = hk and Xsc is the mean-free-path, related by A^c = 7r/2(t'^r). From the 
band dispersion or E versus k relation we define the Fermi wavevector as: 



E'f - EIJA 



hvo 



(10) 
(11) 



where Egap/2+aGqVG marks the position of the Fermi level (Ef) relative to the Dirac-point 
and n5 is the electron density. Relating Eqs. [inland [111 "we can express a voltage and 
bandgap dependent electron density to determine the mobility 

qr qXsc qXsc 



(12) 



m* m*v hkp 

It is thus clear that the mobility depends on the value oi kp- As we vary the bandgap 
of graphitic systems (epi-G, BLG, s-G, or GNR), the variation in kp (equivalently, Ep) 
depends on what parameters are being held constant in the process. To start, let us assume 
the scattering length A is independent of energy, so that we're effectively working in the 
ballistic limit. At this point, we can assume the electron density is constant while the 
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bandgap is being opened, so that kp is constant and the mobihty does not change. However, 
possibly a more suitable metric is the gate overdrive Vg — Vr, which ultimately determines 
the charge density too using Ug = Cg{Vg — Vy), where C^^ = C'^^ + involves both oxide 
and quantum capacitances. In the limit of small density of states {Cq <^ Cox) at smaller 
bandgaps, the gate overdrive is the quantity that is controlled externally, and this changes 
Ug as the quantum capacitance proportional to density of states increases with energy. We 
then get 

aGqVGiaGqVG + Egap) 



ris 



q\ 



(13) 
(14) 



h^y^Tns{aGqVG,Egap) 

where the gate transfer factor is og = Cqx/Cy, and Ce is the equivalent capacitance of 
a three-terminal device including its quantum capacitance. The fundamental mean-free- 
path (Asc) can be approximated semiclassically from the conductance = q^D^Ep)!) ^ 

{Asc^fI with kp defined in Eq. [11] and V being the diffusion constant. Single layer 

h 

graphene (SLG) can have A^c on the order of microns In addition to A^c, a more complete of 
mean-free-path (A) from Eq. [H] would include scattering due to charged impurities, rough- 
ness, and possible phonons from interfacial materials given that phonons native to graphene 
are inherently suppressed, (1/A = l/\sc + ^/ impurities + ^/Kough + 1/Ap/x). 
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FIG. 10. Decrease in mobility (for fixed gate overdrive) between various graphitic materials as well 
as various electron densities 



For a fixed gate overdrive, the mobility even for a ballistic device decreases with bandgap, 
primarily due to the asymptotic constraint that pins the band structure to a high energy 
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linear dispersion. We emphasize that this trade-off arises independent of any reduction in 
scattering length Xsc through the bandgap opening process. The low effective mass of graphene 
arose from its sharp conical bandstructure in the first place, so that opening a bandgap 
without removing the higher energy conical dispersion invariably makes the carriers heavier. 

With respect to digital switching applications, the importance of the above trade-off 
cannot be overstated. The mobility ultimately determines the switching speed through the 
ON current, while the bandgap relates to the ON-OFF ratio. For cascaded devices, it is 
also worth emphasizing that the ON-OFF ratio needs to be computed at high bias, as the 
drain and gate terminals in regular CMOS like cascaded geometries are connected to the 
same supply voltage. Increasing the ON-OFF ratio by increasing the bandgap is predicted 
thereby to reduce the switching speed. We must therefore evaluate GNRs on this entire 
fi — Eg curve rather than at an isolated point on this 2D plot. Eq. [H] elegantly relates Egap 
and /i for various A's. Analyzing the three parameters Egap — fi — X simultaneously allows 
us to project the performance of graphene derivatives and compare against other common 
semiconductors as seen in Fig. [TOl 

18.3. THE THIRD LEVEL ACTIVE THREE-TERMINAL ELECTRONICS 

We now move outward to the next circle on the Y-chart, towards three terminal active 
electronic devices. Our main focus will be on a class of patterned device-interconnect hybrids, 
where we see certain notable advantages mainly on the electrostatics and the contact barriers, 
but challenges with the small bandgap show up as band-to-band tunneling and modest ON- 
OFF ratio. 

18.3.1. Wide Narrow Wide (WNW) - All graphene devices 

The analyses from the previous sections set the platform for evaluating the I-V charac- 
teristics of GNR devices in presence of a third gate terminal. Our lesson from section 18.1 
indicated that the GNR metallicity is primarily set by its ribbon width, showing that one 
might be able to monolithically pattern a wide-narrow-wide all graphene device that flows 
seamlessly from metallic channels to semiconducting interconnects. Experiments have in 
fact shown the ability to carve out GNRs using either chemistry or nanoparticle mobilities 
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that snip the sheets almost perfectly along their C-C bonds. GNRs as thin as 1 nm with 
perfect edges have been manufactured chemically [2^. It is thus interesting to query what 
the device level advantages of such a monolithically patterned GNR would be. We will later 
discuss the circuit level ramifications. 
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FIG. 11. WNW dual gated all graphene device, showing local E — ks (top), top view (center) and 
side view (bottom) with the device parameters listed. 

The structure of an imagined WNW graphene nanoribbon field-effect transistor (GNR- 
FET) is shown in Fig. [TTJ The wide regions are metallic and the narrow ones semiconducting. 
There are planar gates both at the top and the bottom, the top ones for gating and the 
bottom ones for electrostatic 'doping' (see figure later for inverters). Let us first discuss how 
we simulate the I-V of one of these WNW devices. 



18.3.2. Solving quantum transport and electrostatic equations 



The calculations we will show couples a suitable bandstructure/density of states for 
the graphene channel with full 3D Poisson's equation for the electrostatics and the Non- 
Equilibrium Green's Function (NEGF) formulation for quantum transport 15|. The 
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wider contact regions are captured recursively by computing their surface Green's functions 
gi^2{E). The corresponding energy-dependent self-energy matrices Y^i^2{E) = Ti_2 5'i,2Ti2 
project the contact states onto the channel subspace, where the r matrices capture the 
bonding between the contact and channel regions. In order to capture the interfacial chem- 
istry properly, we extend the device a couple of layers into the wider regions and calculated its 
Hamiltonian H matrix. The Coulomb matrix U is computed using the method-of-moments, 



described below 



2l|. 



From the above matrices, the retarded Green's function G = {ES ~ H — U — T,i — S2)~^ 
is computed, and thence the charge density matrix p = J dE GT,^"^ 0^271, whose trace gives 
us the total charge. S*" = (Fi/i + r2/2) in the simple limit where the only scattering arises 
at the contact channel interface. In the previous equation, ri^2 = "^(^1,2 — ^12) ^^e 
contact broadenings (the matrix analogues of the injection rates 71 2 introduced in section 
18.2.1), while fi^2{E) = ^ e^^"^'^-^"'^''^'^] represent the contact Fermi-Dirac distributions, 
with 2 being the bias-separated electrochemical potentials or quasi-Fermi energies in the 



contacts 



161 ]. The charge density matrix is then used to recompute the Coulomb matrix 



U self-consistently through Poisson's equation. Finally, the converged Green's function is 
used to compute the current / = {2q/h) J dE T{E)[fi{E) — f2{E)], where the transmission 
T{E) = trace{TiGT2G^) [jj. 

Let us now get into a few details on the 3D Poisson equation we solve, using the method 
of moments (MOM) numerically. MOM captures the channel potential by setting up grid 
points on the individual device atoms with a specific charge density Sud, and on the contact 



atoms with a specific applied voltage (pc 
for Device, we get 



2l|. Using the notations 'C for Contact and 'D' 




Laplace Single Electron Charging Energy 



where we imply vector notations for the potentials (j) and matrix notations for the Coulomb 
kernels U. And is calculated relative to its neutrality value Nq by tracing over p above, 
while Nq is calculated analogously, while grounding all the contact potentials (this would 
depend on the workfunction of the contacts, as in MOS electrostatics). The matrix elements 
in U need to be computed with the correct dielectric constants. Let us describe it in the 
simpler case with a dielectric constant k for the top gate and a dielectric constant unity for 
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the bottom (trivially generalized to multiple dielectrics). Using the method of images, 



f/(fi,f2) 



4vreoei 



£2 - ei 



1 



|fi - 



(in the same medium) 



27reo(ei + €2) In - f2| 
where is the image of the charge at r2 



e2 + ei, . 

in different media) (16) 
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23| . A tricky point is to avoid the infinities 
X2 and Hi = 2/2 • We can avoid these 



at the onsite locations, for instance, when xi 
using the Mataga-Nishimoto approximation, where we replace terms like l/|ri — "^2! with an 
atomistic correction 1/ a/Iti — r2P + a^, with the cut off parameter a adjusted to represent 
the correct onsite Coulomb (Hubbard) charging en ergy given by the difference between the 
atomic ionization energy and the electron affinity 24 1. 

Let us now discuss the observed electrostatic characteristics in the WNW device, which 
explains the geometric advantages of this particular structure. 



18.3.3. Improved electrostatics in 2-D 




FIG. 12. (Left) The two-terminal potential shows the vanishing fields near the channel, implying 
the superior gate control and the improved short-channel effects with the 2D contacts. (Right) 
The 3D potential shows the non-linear flat potential in the middle of the channel. 

We simulate a device patterned monolithically from a two-dimensional sheet of graphene 
with a wide dilution of widths from the source and drain contacts to the active channel 
region. Simulated WNW (35-7-35) GNRFETs compose of (7,0) armchair graphene nanorib- 
bon (GNR) narrow regions for the channel and (35,0) armchair GNR regions for the contact 

20 



and interconnect regions. A metallic gate is placed on top of the channel region, while a wide 
grounded substrate is placed at the bottom of the channel. For calibration with the conven- 
tional CMOS technologies, the unique two-dimensional (2D) contacts of the GNRFET are 
replaced with 3D bulk metal contacts (whose surfaces act as parallel plate capacitors) for 
the same device, gate, and dielectric geometry. 

A particular advantage of the WNW structure is the low capacitance of the 2D source 



drain contacts [25|. In a conventional MOSFET, the gate electrode needs to compete elec- 
trostatically with the source and drain for control of the channel charge. Indeed, a majority 
of developments in transistor technology over the last few decades have concentrated on 
making the field lines gate controlled rather than source/drain controlled. This is becom- 
ing harder with aggressive size scaling. The 2D source and drain contacts with a top gate 
makes the S/D capacitances lower, as they can only influence the channel through their 
fringing fields. Note that a 2D side gate geometry, as advocated in many device designs, 
would eliminate that electrostatic advantage, as the gate needs to compete with the S/D 
electrodes. 

As the channel length gets shorter with the aggressively scaled technologies, the 3D 
contacts start to influence the channel potential as their surfaces act as parallel capacitor 
plates flanked by the insulator at the top and bottom. In the case of the monolithically 
patterned 2D GNR contacts, the charges on the contact surfaces are line charges so that the 
applied source-drain field decays into the channel, creating a non-linear channel potential 
even in the absence of a gate (Fig. [121 left). Moving on to a three-terminal, dually gated 
structure. Fig. [12] (right), shows that the gate contact holds the channel potential flat against 
the action of the drain, thereby reducing short-channel effects. 

Figs. [T3| [14] show that for the same channel geometry, the top gate with 2D side contacts 
has the largest capacitance, followed by the top gate with 3D side contacts and finally 
the lowest gate to drain capacitance ratio is obtained when all electrodes are co-planar. 
The corresponding field line diagrams are also shown in these figures. Note also that in 
addition to the source, drain and dual gate electrodes, one needs to worry about the quantum 
capacitances, which are automatically included from our density matrix calculations that 
enter Poisson's equation. 

The capacitance ratio can be extracted by plotting the channel transmission (T) for two 
scenarios: maintaining a constant drain voltage (V^) while sweeping gate voltage (V^) , and 
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3D Source-Drain 



2D Source-Drain 





FIG. 13. Comparison of planar source drain vs 3D source-drain. Denser field lines on the channel 
from the 3D contacts correlate to stronger source coupling and DIBL. For the given material and 
geometrical parameters listed in Fig. [TT\ the Cg/Co ratios are 4.95 and 5.80 respectively. Top 
and bottom gates were grounded while the source was simulated with a potential of 0.3V and 
conducting channel had a potential of O.IV 



analogously, maintaining a constant Vg while sweeping V^. Sweeping the Vg creates a larger 
energy shift in the transmission of the GNRFET channel than the sweeping of the V^. From 
the shifting rates of these transmissions and the charge density calculations from the MOM, 
we can extract the capacitance values of the contacts. With shifts in transmission plots, 
we once again find that 2D contacts indeed help the gate exercise superior control over the 



channel 
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We will now explore the effect of the improved short channel effect on the computed I-V 
characteristics. 



22 




Channel 



FIG. 14. Comparison of top vs side gate. Denser field lines from the top gate ensure better 
gate control which is reflected by larger gate capacitance. For the given material and geometrical 
parameters listed in Fig. [TTl the Cg/Cd ratios are 5.80 and 3.82 respectively. Gates were biased 
at 0.4V , while the conducting channel had a potential of O.IV. 

18.3.4. Three terminal I-Vs 

The computed three terminal I-Vs (Fig. [T5l) show excellent short channel effects, at least 
over a small voltage range given by the bandgap. Plotted vs gate voltage, the current shows 
excellent saturation characteristics with a large output impedance. Plotted vs gate voltage, 
the current shows little drain bias dependence (so-called drain induced barrier lowering or 
DIBL). Taken together, the curves signify that the device electrostatics in the geometry is 
nearly ideal, making the outputs relatively robust with process variations. It is interesting 
to note that instead of enhancing the gate capacitance as in regular CMOS devices, the trick 
in WNW devices has been to reduce the source and drain capacitances in comparison. 

The simulations results of the model in the Fig. [T5] demonstrate a Subthreshold Swing 
(SS) of 84.3 mV/dec and a Drain-Induced Barrier Lowering (DIBL) of 24 mV/V. We note 
that unless otherwise specified all simulations refer to material and geometrical parameters 
shown in Fig. [11] The value of the DIBL and the SS can be further improved by increasing 
the length of the channel (currently 1:8.6 ratio of Hf02 thickness to channel length). These 
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values calculated are better (smaller) than the estimated values of DIBL = 122 mV/V and 
the SS = 90 mV/dec for the double gate, 10 nm scaled Si MOSFETs |26|. Also in addition 
to showing improved short-channel effects, the GNRFET structure with the 2D contacts also 
shows controlled switching behavior. The on-current (Ion) of the system equals to 2670.62 
Aj fim with the off-current (loff ) set at 4.07 A/ fim; thus giving a lon/Ioff ratio of 656. The 
ON-OFF ratio, however, ends up being modest, and is a critical challenge in GNRFETs, 
especially in the light of its seemingly inverse relation with the charge mobility (section 
18.2.2). 
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FIG. 15. I-V curves for a n-type GNRFET confined to create a large bandgap (in this case, a 
(7,0) armchair GNR with a bandgap nearly 1 eV). Such an extreme geometry postpones the onset 
of band-to-band tunneling. More importantly, the point of the I-V is to show the effect of better 
electrostatics which is independent of bandgap issues - resulting in a high current saturation, low 
DIBL and SS. 

With the scaling of the channel length, the short-channel effects started to have a huge 
influence on the device parameters such as the DIBL and SS. As the channel lengths get 
shorter, the DIBL and SS of the device increases due to the Cg/Cd ratio of decreasing with 
length. The line charges with the 2D contacts endow the gate with more control over the 
channel and interface states compared to the 3D contacts by lowering the drain capacitance. 



18.3.5. Pinning vs. Quasi-Ohmic contacts 



In today's semiconductors, Ohmic contacts are a desired to help achieve linear and asym- 
metric I-V characteristics. The potential profile inside the channel can be influenced by 
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increasing the drain-source voltage (Vds) or the gate voltage (Vgs). For carbon nanotubes, 
this has been a particular challenge, as the metal carbon bonds at the ends have predom- 
inantly created Schottky barriers 27|. In our WNW geometries, since the bulk metal 
contacts are relegated to the ends of the device array, the bonding configuration near the 
wide-narrow interfaces are controlled by C-C covalent bonding. As our simulations show, 
this seems to promote a quasi-Ohmic behavior. The better bonding increases the decay 
lengths of the corresponding metal-induced gap states (MIGS) entering from the wide re- 
gions. The partial delocalization reduces the single-electron charging energy (that enters 
through our MOM treatment), thus making it harder for the contact regions to pin the 
Fermi energy and reducing the effectiveness of the Schottky barrier. 

Schottky barrier FETs behave qualitatively different from MOSFETs In the latter, an 
applied gate bias reduces the channel potential and controls the thermionic emission over 
the volt age- dependent interfacial barrier. In the former, the gate reduces the thickness of 
the Schottky barrier and controls the tunneling of electrons through a voltage-independent, 
pinned barrier height. The question is what the potential profile looks like in the channel, 
and whether the contact MIGS are effective in pinning this potential adequately. 

As seen in the Fig. [161 the lowering of the potential throughout the entire graphene 
channel region with applied gate bias is a characteristic of the regular ohmic contact FETs 
rather than the Schottky barrier FETs, whose potentials would otherwise be pinned to the 



28|. 



midgap by the charging of the interfacial states 

The MIGS due to the tail ends of the metal states in the contacts, leak in the semi- 
conductor. Even with this 2D contact geometry, the MIGS will be present because of the 
contact-channel interfaces 25|. Our WNW all-graphene structure can filter these quickly 



decaying states, resulting in no significant contribution to the electron transmission. In the 
case of our device with the channel length of 8.66 nm, the MIGS do not travel all the way 
from source to the drain, but only extend approximately 0.7 nm into the semiconducting 
channel (Note a typo in one of our earlier papers, where we wrongly quoted this as 0.07 nm) 



25| . The decay length of these MIGS can be calculated by plotting the wavefunction of 



the channel electrons at specific energies, as well as by evaluating the complex E-k diagram. 
The intensity of these MIGS at a given distance x can be calculated by using the equation 
lQ*e^~^^'^'^\ where Iq is the intensity of MIGS at the interface and A is the decay length. 
Note that issues similar to those discussed here have been discussed in the context of 
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FIG. 16. At Viis = O.OV and Vds = l.OV, variation of channel potential with different gate voltages 
shows no barrier pinning at the contacts, implying Ohmic contacts. 



pentacene molecules with CNT contacts. While CNTs would offer even better ID electro- 
static gains, a trade-off arises with the increasing series resistance in CNTs due to a paucity 
of modes. For GNR source/drain analogously, we will need to imagine wide blocks simul- 
taneously contacting many GNR devices, so that the contact resistance is minimized by 
extending its width. 

We now have all the tools to compute three-terminal I-Vs in graphitic structures, doing 
full justice to the complex electrostatics. Let us now see how this influences the circuit level 
performance metrics of GNRs. 



18.4. THE PENULTIMATE CIRCLE: GNR CIRCUITS 

Moving to the next level of abstraction and the last rung of our Gajski-Kuhn Y-chart, 
we focus on circuit level design issues when integrating GNR interconnects or devices. We 
start by discussing potential circuit topologies where we advocate a truly all graphene cir- 
cuit, which is a logical extension to our WNW all graphene FET analogue. Next, instead 
of delving into an array of circuits we present fundamental circuit response design issues 
revolving around some basic GNR based logic element, namely the inverter. The inverter is 
the most basic logic circuit with the sole purpose of converting a logical 1(0N) to a logical 
O(OFF) and vice versa for opposing input voltages. In the process of characterizing a GNR 
inverter we touch on interconnects and finally discuss drivability between connected GNR 
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logic elements. 




FIG. 17. The various building blocks for all graphene circuits 



18.4.1. Geometry of an all graphene circuit 

An attractive feature of graphene for device engineers is its planarity, which stems from 
its sp^ chemical bond hybridization. Graphene's atomic flatness is compatible with existing 
lithographic device fabrication techniques established for CMOS. Furthermore, exploiting 
how chirality in graphene influences electronic properties we present an array of building 
blocks for all graphene circuits (Fig. [T7|) . In the previous circle on our Y-chart we explored 
an all armchair GNRFET. 

Circuit level enhancements start at the device level. Two important device performance 
metrics for digital circuit performance are ON-OFF current ratio and intrinsic gate propaga- 
tion delay. GNRFET ON current scales proportionally with width, while the OFF current 
goes as e^^^'''^ or g'^/'^^^^ where W is the width. To achieve manageable OFF currents for 
digital applications, GNR widths must be scaled within the sub-lOnm regime to avoid in- 
creased in static power dissipation and poor ON-OFF ratio. However GNR scaling has little 
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influence on propagation delay. 

Propagation delay is defined as CjNVdd/ Ion, where Cjn is the intrinsic capacitance, 
Vdd is the supply voltage and Iqn is the saturating ON current at Vgs=Vdd- Increase in 
channel width increases Iqn as threshold voltage decreases, which simultaneously drives 
up the OFF current exponentially. However Cjn = {I /Cox + ^/Cq)^^, which also scales 
with width balances any improvement in propagation delay from scaling Jo at in graphene. 
From the perspective of a single GNR device there are not many options to improving the 
propagation delay. However a more useful context to address this issue would be on the 
circuit level where we have multiple logic elements. 

3D metal contact Devices in Series 3D metal contact 




metal contacts 



FIG. 18. Geometry of all graphene circuits 

On a circuit level the complexity of propagation delay is in the load capacitance, which 
is now the sum of all capacitances at the output including the combinational drain capaci- 
tance from the previous inverter, interconnect capacitance(Cj„fer'connect); and next load gate 
capacitance. Using the building blocks in Fig. [T7] we present an all graphene circuit shown 
in Fig. [TBI with cascaded GNRs in parallel on a semiconducting substrate with separate 
split gates for electrostatic doping the regions n and p-type. The cascaded GNRs could be 
separated by a high-k dielectric or even boron nitride in hexagonal lattice, which has the ad- 
vantage of being atomically fiat and the absence of dangling bonds 17| makes it less likely 
to carry adsorbents that could degrade the device. The advantage here is that ON-OFF 
ratio is held constant, while the increased ON current and capacitance which scales propor- 
tionally with the number of cascaded GNRs dilutes the parasitic interconnect capacitance, 
improving propagation delay and circuit performance. 

Good quality graphene sheets have been made viable by current advancements in wafer- 
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scale and pattern transfer techniques 18|, |29|]. However full realization of an all graphene 



circuit with various GNR interconnects and devices rely on the ability to pattern GNRs to 
narrow enough widths to produce a sizable bandgaps. Planar lithographic techniques are 
prone to edge roughness, while various chemical methods have had the most success in cre- 
ating chemically precise GNR edges, but their applicability to scalable device level processes 
remain to be seen Q jsO] 31 1 32 1. While roughness helps to make GNRs insensitive to 
chirality, we need further simulations to see how atomic fluctuations in the widths influence 
the corresponding threshold voltages and ON/OFF currents, an issue critical for the overall 
reliability of GNR circuits. 



18.4.2. Compact Model Equations 

To simulate the performance of such a circuit, let us first outline a compact model. 
This will require us to outline (a) an equation for the bandstructure that includes effects 
due to edge strain and roughness, (b) an equation for the scattering length that depends 
on the phonon spectrum and edge roughness, (c) equations for the 2D electrostatics due 
to the source and drain contacts, and (d) the resulting I-Vs obtained by integrating the 
transmission over the relevant energy window. 

To recap, the bandstructure of GNRs, including edge strain, can be written in a tight- 



binding form a.s E = ±^Eljy + h^vlkP'. Specific expressions for Ecy and Vq for variously 
strained graphitic materials exist in Ref. j9|. 

The next term is the scattering A^c, which is related to the scattering time through an 
angle averaged geometrical factor and the overall Fermi velocity. The scattering time is 
extracted from Fermi's Golden Rule. For short range scattering by edge roughness and 
phonons, r^c oc while for long ranged unscreened Coulomb scattering, r^c oc \E\. 



Explicit expressions exist in the literature 33| 34 1. 



The tricky part that does not exist in the literature are the electrostatic capacitances 
associated with the 2D electrostatics from the planar source and drain contacts, competing 
with the top and bottom gates through their individual dielectrics. We are in the process 
of extracting formulae based on knowledge of planar micro-strip line electrostatics, with 
geometrical factors calibrated with our numerical MOM solutions for a variety of geometries 
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Once we have the electrostatic, band and scattering parameters, we can then use Eq. [T] 
to extract the I- Vs. For energy-independent A^c, this was aheady shown earher. We will 
generalize it to various scattering configurations in our future work. 

We thus have a comprehensive compact model that captures the chemistry and bandstruc- 
ture, scattering, electrostatic and transport parameters needed for our circuit simulations. 
We will report one example here, and report further results in our subsequent publications. 



18.4.3. Digital circuits: 

Static complementary CMOS gates utilize pull-up (PUN) and pull-down (PDN) networks 
to achieve low power dissipation and large noise margin in logic circuits such as the inverter, 
NAND, and NOR gates. CMOS logic circuits are composed of some series and parallel 
combinations of n and p-type FETs. An inverter illustrated in Fig. [19] is the simplest logic 
element and the focus of this section of the review. 

Inverter Characteristics 




FIG. 19. GNR inverter geometry and voltage transfer curve. This inverter design uses the WNW 
(metal-semiconductor-metal) all graphene structure for pull-up and pull-down networks. In this 
design CNT interconnects make direct contact with device level graphene. CNT/graphene interface 



has been experimentally demonstrated by Fujitsu Laboratories Ltd 
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When the input into the common gate is Vin=0 the p-type FET (PUN) is active while 
the n-type FET (PDN) is cut-off, hence the circuit will pull the output voltage up toward 
the supply voltage (Vdd) or high, Vout=l. Likewise when Vin=l, n-type FET is active 
and p-type FET is cut-off pulling the circuit down toward ground, Vout=0. Usually it 
is impossible to pull-up or pull-down to exact values of 1 or 0, so threshold voltage and 
tolerance are designed for each circuit to help distinguishing between these two logic levels. 
Circuit designers allow some tolerance in the voltage levels used to avoid conditions that 
generate intermediate levels that are undefined. For example, to 0.2V on the output can 
represent logic (0) and 0.3 to 0.5V can show (1), making the 0.2 to 0.3V range invalid, not 
metastable, since the circuits cannot instantly change voltage levels. 

The voltage-transfer curve (VTC) of an inverter circuit captures the DC or steady-state 
of specific input versus output voltages and provides a figure of merit for the static behavior 
of the inverter. VTCs for logic circuits provide information on operating logic-levels at the 
output, noise margins, and gain. Ideally we want the VTC to appear as an inverted step- 
function, indicating precise switching between the on and off stages, but in real devices there 
is a continuous transition between on and off states. From the VTC we can extract noise 
margin (Fig. [20l) . which provides a measure of circuit reliability and predictability. Biasing 
outside the noise margin puts the logic circuit in an unpredictable state. Circuit designers 
want to maximize the noise margins. 




FIG. 20. showing the importance of balancing CMOS transistor sizes to achieve equal high and 
low noise margins (NM). The noise margin is graphically represented by the largest square that fits 
inside the enclosed space outlined by normal and rotated VTCs. 



A significant advantage of graphene is its intrinsic electron-hole effective mass symmetry. 
In the absence of extrinsic doping a graphene based FET the I-V characteristics for n and 
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p-type conduction would be the symmetric. However, asymmetry can be introduced into the 



system through charge-transfer doping 



17| (e.g. Fig. [8]) or by contact induced doping [3 



Significant screening of charge impurities in the substrate should bring Fermi level closer to 
its intrinsic value at the Dirac or A'— point, therefore recovering symmetric n and p-type I-V 
characteristics. On a circuit level this symmetry means the response of PUN and PDN would 
he equal and opposite, which is important for circuit reliability, and not to mention ease of 
circuit design. In conventional Si-CMOS logic circuits, the asymmetry in the electron-hole 
effective mass is compensated by scaling the physical width of the p-type FETs in the PUN 
so the I-Vs are equal and opposite with the PDN. Graphene's natural electron-hole symmetry 
would allow circuit designers to bypass this design issue. 

A major impediment to GNR based logic circuits is its narrow bandgap ( ^ 200meV), 
as the device elements in the PUN and PDN are prone to sub-threshold leakage from band- 
to-band tunneling. The two-fold effect on an GNRFET-based inverter where the channel 
has a narrow bandgap is demonstrated in Fig. [21] The first effect is a large voltage swing of 
approximately 0.4V. The second effect is a significantly diminished noise margin. Band-to- 
band tunneling in narrow bandgap GNRFETs prevents either the PUN or the PDN from 
completely cutting off when its complement network is active. 




FIG. 21. Comparison of VTC curves for narrow bandgap GNR and 45nm CMOS Technology. 
Narrow bandgap GNRFETs will be more susceptible to noise than CMOS due to smaller noise 
margins. 

Fig. HH] shows the physical layout of a functional graphene inverter composed of WNW 
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P-type and N-type GNR device arrays and the voltage transfer curve. The inverter voltage- 
transfer curve and gain can be calculated readily from the current-voltage characteristics. 
As expected the gain of the device determined by the electrostatics, geometrical parameters, 
and mobilities which ultimately determine the P and N-type GNR transconductances. The 



VTC above with gain of 4 is derived from the I-V s. 
by using the methodology described in detail in 138 



lown in Fig. for the 8.66nm device 
These I-Vs generated in SPICE can 



be used to simulate other complex layouts such as NAND or NOR gates shown in Fig 
(The results of these logic gates will be reported in future publications). 

Propagation delay can be measured by pulsing the input voltage between and 1 and 
observing the output transient response. The transit time for a GNRFET is approximately 
L/v, where L is the length of the channel and v is an energy dependent velocity defined in 
Eq. [3l Intrinsic and extrinsic device level scattering mechanics could also influence transit 
time. However, a cascade of inverters or some other logic elements in series, the load ca- 
pacitance between each logic stage typically dominated by Cmterconnect would be responsible 
for the majority of the delay. The GNR circuit layout we presented earlier and show in 
Fig. HH] takes advantage of non-scaling parasitic interconnect capacitances, lowering delay 
and improving performance. 

Beyond individual logical elements (ie., inverter, NAND, NOR), an important CMOS 
circuit design parameter is fan-out, which estimates the number of logic stages or CMOS 
gates that can be consecutively driven before signal attenuation is no longer tolerable. Past 
the maximum fan-out a repeater or amplifier is necessary to drive subsequent logic stages 
in a circuit. The maximum fan-out scales proportionally with propagation delay; therefore 
circuits designed for low frequency applications have a larger maximum fan-out compared to 
circuits designed for higher frequency applications. If graphene is to indeed follow the MOS- 
FET and CMOS paradigm fan-out would be important circuit design trade-off to consider 
and a topic we will discuss in an upcoming work. 



18.4.4. Physical domain issues: monolithic device-interconnect structures 

The bane of many CNT-based circuit ideas is the degradation due to the dominant 
aspect of the Schottky contacts between the devices and interconnect. Luckily, GNR circuits 
can avoid this problem by using the same sheet of graphene for both active devices and 
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interconnect, as seen in Fig. [22] 

Thus, for such monohthic GNR device/interconnect structures with Ohmic contacts, the 
behavior can be close to the ideal device predicted by simulations. Of course, contacts 
to metal are still required, but at a coarser granularity, thus the expected degradation in 
performance will be proportionally smaller. Contacts to metal are necessary for several 
reasons; first, the FET gates and drain/source electrodes cannot share the same graphene 
sheet, so connecting the output (drain) of a device to the input (gate) of the next will require 
a contact, second, topological requirements for connecting complex circuits can rarely be 
mapped to a planar graph, and non-planar graphs implicitly require more than one layer 
of interconnect, thus contacts, and for a really large circuit a large number of interconnect 
layers (~ 10 for a modern CMOS circuit). 



18.5. CONCLUSIONS 

The future of microelectronics relies on continually scaling the critical dimensions of bulk 
CMOS technologies. The semiconductor industry faces serious challenges in this respect, 
due to a host of technical and economic constraints. One way to mitigate this is to use 
novel architectures, combining inherently scalable top-down techniques, such as lithography, 
with novel bottom-up fabrication approaches, such as self-assembly. An alternate way is to 
look for novel channel materials beyond silicon - a strong candidate being another group IV 
element, carbon, graphene being one of its distinguished allotropes. 

Until recently, much attention has been focused on CNTs. CNTs have already demon- 
strated excellent intrinsic performance, high gain, high carrier mobility, high reliability, and 



are almost ideal devices in themselves 



39| 40| j41|. Unfortunately, there seem to be no 



practical, scalable solutions for arranging multiple CNTs on a substrate uniformly with a 
small pitch, needed to deliver adequate current for fast switching 3l| 42| . Neither are there 
clear approaches for contacting CNTs to interconnects or to each other at low impedance 
to realize complex circuits using large scale fabrication schemes. In fact, CNTs frequently 



encounter Schott 
their reliability 



iv barriers at the contacts that control the tunneling electrons, hampering 
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In contrast, graphene's planar profile makes it amenable to well-established planar fab- 
rication techniques for silicon CMOS devices. Its mobility can reach up to two orders of 
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FIG. 22. GNR NAND layouts with electrostatic doped back-gates and interconnecting 'vias' be- 
tween multiple levels. CNT/graphene interface has been experimentally demonstrated by Fujitsu 



Laboratories Ltd 
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magnitude above silicon, and the ability to engineer its bandgap with width alone points 
to the feasibility of all-graphene devices that can exploit covalent bonding chemistry at the 
contact and better inherent electrostatics to allow more traditional, MOSFET like gate con- 
trol mechanisms. However its bane seems to be its metallicity. As we saw earlier, graphene 
is naturally a zero-band gap semi-metal, and attempts to open a bandgap, such as through 
strain, quantization or field asymmetries have limited the bandgaps to < 200 mV [20[. In 
a regular field effect transistor application, such a bandgap translates to an ON-OFF ratio 
of ~ 70, inadequate for digital logic 4^. This seems consistent with experiments on high 
speed graphene transistors, which show I-Vs that are essentially quasi-linear. There is a lot 
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of activity also on opening bandgaps in graphene, although it seems that this may reduce 
the mobility that made graphene a promising electronic device in the first place. 

Despite the significant challenges, graphene's high carrier mobility and fast switching 
speeds make it widely studied. This includes revolutionary applications such as based on 
charge focusing, graphene spintronics and excitonic condensation of pseudospin states 45 |. 
At the same time, there is wide activity on graphene based conventional electronic devices, 
specifically, RF devices and CMOS switches. It is not clear what the prospects of GNRs are 
vis-a-vis switching, given its low bandgap and its seemingly fundamental mobility-bandgap 
tradeoff. However, as we have argued here, there still are a few notable advantages to 
using graphene geometries, namely, (i) a natural electron-hole symmetry that helps with 
inverter design; (ii) its convenient placement between ID and 3D, so that it offers distinct 
electrostatic advantages without picking up too much series resistance arising from a paucity 
of modes; (iii) the tunahility of its electronic properties primarily by controlling its width. 
Experiments are still emerging on these fronts, and it remains to be seen whether one can 
put these advantages to good use. 

In this chapter, we made a first attempt to transverse the Gajski-Kuhn Y chart adapted 
to graphene, spanning physical, structural and behavioral domains for graphene. We ended 
with monolithic all graphene circuits and device/interconnect geometries, ultimately stop- 
ping short of the outermost circle that covers a system level topology and behavior. A 
key tool developed along the way was a compact model for GNRs that capture chemical, 
geometric, electronic, electrostatic and transport properties through a set of simple parame- 
terized equations calibrated against a few emerging experiments. Such a Y-chart is the first 
step towards a holistic approach that we hope will catapult graphene from the domain of 
fascinating physics and chemistry to technologically relevant electronic applications. 
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